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PREFACE 


This  semiannual  technical  report  describes  the  work,  performed 
under  Contract  Number  F  44620-75-C-0088  at  Harvard  University,  during 
the  period  October  1,  1976  -  March  31,  1977. 
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Mr.  P.  Liu,  Graduate  Student 

Material  in  previous  reports  (see,  for  example.  Cruft  Laboratory 
Technical  Report  665,  October  1976),  is  not  duplicated. 

The  views  and  conclusions  contained  in  this  document  are  those  of 
the  authors  and  should  not  be  interpreted  as  necessarily  representing  the 
official  policies,  either  expressed  or  implied,  of  the  Defense  Advanced 
Research  Projects  Agency  or  the  U.S.  Government. 
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A.  REPORT  SUMMARY 

High  power  excimer  lasers  operating  in  the  ultraviolet  region  of 
the  spectrum  have  revealed  limitations  imposed  by  failure  of  window 
materials  and  coatings.  More  reliable,  quantitative  data  on  two  photon 
absorption  cross  sections  and  other  nonlinear  processes  in  ultraviolet 
window  and  coating  materials  are  needed. 

We  have  developed  at  Harvard  University  over  the  past  several 
years,  various  techniques  to  study  quantitatively  breakdown  thresholds 
and  other  nonlinear  phenomena.  Reproducible  single  picosecond  pulses 
with  Gaussian  temporal  and  spatial  profiles  are  available  from  a  Nd-Yag 
system.  We  have  used  pulses  at  the  fourth  harmonic  frequency,  corres¬ 
ponding  to  the  ultraviolet  wave  length  X  =  266  am,  to  measure  the  two- 
photon  absorption  cross  section  in  a  number  of  alkali  halide  crystals. 

The  results  are  detailed  in  Section  B  of  this  report.  The  intrinsic  two- 
photon  absorption  coefficients  lie  in  the  range  1  -  4X10*^  cm/Mega  watt. 
Details  may  be  found  in  Table  1  of  Section  B.  The  cross  sections  appear 
to  be  about  one  order  of  magnitude  smaller  than  previously  reported  for 
KC1  and  RbBr.  Lithium  fluoride  with  a  bandgap  larger  than  2  Au,  does 
not  show  any  detectable  two-photon  absorption  at  266  am.  The  data  will 
be  extended  to  other  wavelengths,  notably  the  third  harmonic  at  355  nm, 
and  to  other  oxide  and  fluoride  window  and  coating  materials. 

A  ruby  laser  system,  coupled  with  a  dye  laser  pumped  by  the 
second  harmonic  pulse  of  the  ruby  frequency,  has  been  used  to  measure 
the  two-photon  absorption  cross  section  in  TiO^,  GaP,  CdS  and  SrTiO 
A  novel  method  has  been  developed,  in  which  the  two  photon  absorption 
cross  section  is  compared  with  known  attenuation  due  to  the  inverse 
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raman  effect  in  cyclohexane.  Details  are  described  in  Section  C  of 
this  report,  which  is  essentially  a  manuscript  submitted  for  publication 
in  the  Physical  Review.  This  work  will  be  extended  to  higher  frequencies 
and  materials  with  wider  bandgaps.  It  would  be  of  particular  interest  to 
vary  the  dye  laser  frequency,  so  that  the  two  photon  absorption  cross 
section  can  be  scanned  across  the  intrinsic  band  gap. 

Finally,  the  nonlinear  infrared  absorption  in  InSb  is  being  studied. 
The  subnanosecond  CC^-pulsed  laser  facility  permits  the  study  of  non¬ 
linear  effects,  at  intensity,  which  in  longer  pulses  would  invariably 
damage  the  window  material.  The  intrinsic  distortion  of  valence  and 
conduction  band  in  InSb  has  been  studied  theoretically  by  others, 
notably  Tsoar  and  Gersten.  A  nonlinear  behavior  of  the  transmittance 
and  reflectance  of  InSb  during  the  short  pulses  has  been  observed. 

These  data  are  presented  in  Section  D  of  this  report.  Their  interpre¬ 
tation  in  terms  of  solid  state  plasma  theory  is  continuing. 
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B.  TWO-PHOTON  ABSORPTION  OF  ALKALI -HALIDES  AT  266  nm 


Due  to  the  increasing  interest  in  u.  v.  lasers,  it  will  be  useful  to 

know  the  absolute  value  of  absorption  coefficients  of  materials  with  large 

band  gap  like  Alkali -halides  which  have  band  gap  range  from  5.  6  eV  to 

12.  5  eV.  Few  data  were  published  for  Alkali -bromides  and  Alkali- 

iodides  [1-4].  In  this  experiment  the  absorption  of  Alkali -halides  were 

studied  with  the  4-th  harmonic  of  YAG-Nd  laser  at  266  nm. 

The  Alkali -halides  can  absorb  the  266  nm  photon  (4.  66  eV) 

appreciably  through  nonlinear  processes  at  high  intensity.  The  two 

(2) 

photon  transition  rate  W'  is  given  by  second  order  perturbation 
theory  [5-6],  and  Keldysh  model  [7]. 

The  method  used  is  described  elsewhere  [8].  The  attenuation  of 
the  beam  is  given  by 


|L  =  -d-P!2 


(1) 


a  is  the  one -photon  absorption  coefficient,  and  8  is  the  two  photon 
absorption  coefficient,  which  is  related  to  by 


8  =  2 


W 


(2) 


(2) 


The  transmitted  intensity  can  be  found  from  (1) 


I(r,jt,t)  = 


(l-R)2I(r,  L  ,  t)  e_Qr^ 


1  +  8(l-R)I(r,jfc  ,t)(l-e"a^)/<* 


(3) 


where  l  is  the  thickness  of  the  sample.  R  is  the  surface  reflectivity 
the  multiple  reflection  within  the  medium  is  neglected. 
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The  mode-locked  YAG-Nd  laser-amplifier  system  is  described  in 
[8]  and  in  a  previous  joint  technical  report  (ARPA  and  Joint  Services 
Electronic  Program  No.  665,  October  1976).  The  system  gives  a  10  mJ 
single  pulse  at  1.06  pm.  The  beam  profile  was  measured  carefully.  It 
is  Gaussian  both  in  space  and  time.  The  pulse  has  duration  30  ps.  at  full 
width  half  maximum.  The  fundamental  is  frequency  doubled  by  CDA 
crystal,  then  doubled  again  by  ADP  crystal  (Fig.  1).  Both  are  tempera¬ 
ture  tuned  90°  phase  matching  which  gives  no  walk-off.  The  overall 
power  conversion  efficiency  at  100  MW  input  is  20%.  The  beam  profile 
at  266  nm  is  monitored  by  photo-diode  array,  and  fitted  with  a  Gaussian 
curve  with  1/e  radius  0.  68  mm. 

In  order  to  give  absolute  values,  the  beam  is  monitored  carefully. 


A  plane  vacuum  photo-diode  (Fig.  1)  is  used  with  fast  oscilloscope 

to  make  sure  the  beam  is  a  single  pulse.  The  fluctuation  in  the  pulse 
width  is  taken  into  account  by  monitoring  the  energy  of  each  pulse  both 
at  1.  06  pm  and  266  nm  [9].  Consider  the  fundamental  Gaussian  pulse  as 


I(r,  t)  =  Ime 


-(r/d)‘ 


-(t  fr)* 


1/2  2 

The  energy  is  given  by  E^  =  TT  1^  tA,  where  A  =  rtd  .  The  4-th 

1/2  4 

harmonic  has  energy  E,  =  qTT  I  r'A’  where  r'  =  r/2.  A'  =A/4, 

n  m 

and  q  is  the  conversion  efficiency.  We  have 


The  intensity  at  the  sample  is  increased  by  a  pair  of  focusing- 
collimating  lenses.  The  sample  is  put  behind  the  focus  to  avoid  self- 
focusing.  The  focal  length  of  the  lens  is  measured  by  comparing  the 


FIG- 1  BS  .  BEAM  SPUTTERS  D.  PHOTODIODES  F.  FILTERS 
R.  RETICON  PHOTODIODE  ARRAY. 


beam  sizes  before  and  after  the  lenses  system.  With  a  photo-diode 
array,  the  focal  length  is  determined  to  be  24.  7  ±  0.  3  cm. 

In  order  to  determine  0,  the  energy  transmission  coefficient  is 
measured.  For  a  Gaussian  beam: 

Im.-<r/d,2e-<t/T)2(l-R)2e-°'i 

I(r,i,t)  =  - - - — -  (6) 

1  +  01  1-R)  e"(r/d)  e‘(t/T)  (l-e~ai)/a 


T  = 


,out  _  2tt/q®  rdr Im(r,Z,t) 


T  _  ,2  1/2 
I  TTd  TT  T 
m 


/v  -<*  A 

T  ®  e 

2^77^ 

m 


l-e-ai  •'o 


I  u(1+l 


I  (l-R)(l-e 
m 


At  low  intensity,  we  can  get  an  approximate  form  for  the  integration. 


We  find: 


lim  T  -  e~a*(l-R)2 

I  ->° 
m 


ii»  m-  -  i 

t  .  n  d  I  a 

I  -*  0  m 
m 


®  2V2  (1-R)  e 


The  approximation  is  good  for  (0/or)  I  (l-R)(l-e”a,e)  <  1,  which  is 

m 

true  in  the  experiment  for  I  as  large  as  600  MW/cm2.  Also,  the 

m 

one-photon  absorption  coefficient  or  is  small  in  Alkali-halides.  We  get 


lim  T  “  (1-R) 

I  0 


pi 

2  -Jz  ( 1-R) 


(11) 


lim 


m 


d(l/T) 
d  I 

m 


Now,  it  is  clear  that  the  slope  in  1/T  verses  I  .  graph  will  give  the 
P  value. 

I _  is  measured  with  D,  which  is  calibrated  with  respect  to  a 

m3 

thermocouple  with  calibration  accuracy  -  5%. 

Typical  data  is  shown  in  Fig.  2.  The  rather  large  absorption  at  low 
intensity  can  be  ascribed  to  surface  polishing.  From  the  approximate 
value  of  P  found  by  Equation  (11),  a  set  of  curves  are  generated  by 
Equation  (7)  with  different  values  of  P  near  that  approximate  value. 

The  P  is  then  chosen  by  the  best  fit. 

In  order  to  make  sure  there  is  no  artifact,  the  transmission 
coefficients  are  measured  without  sample  and  with  sample  having  band 
gaps  larger  than  2  ftu.  Typical  data  for  LiF  are  shown  in  Fig.  3. 

The  error  involved  in  the  determination  of  P  is  found  to  be  8% 
from  beam  size,  12%  from  the  approximation  of  constant  conversion 
efficiency  g  which  is  intensity  dependent  [10j,  5%  from  calibration  of 

photodiode,  and  the  scattering  of  data  points  which  is  15%  to  30%. 

Thus,  the  uncertainty  is  40%  -  55%. 

Fluorescence  is  observed  in  the  medium  along  the  path  of  the  beam 
[11,  12].  Visible  color  can  be  seen  after  a  few  shots.  In  KC1  the  absorp¬ 
tion  of  5  32  nm  by  the  color  track  is  observed  which  can  be  explained  by 
the  position  of  absorption  band  of  F-center  [13]. 

The  physical  process  involved  seems  to  be  electron  hole  pairs 
generated  by  two-photon  excitation.  They  form  excitons  and  give  out 
fluorescence.  The  exciton  energy  is  given  up  to  form  Cl^  and 
center  [14]. 


F- 


FIGURE 
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By  measuring  the  transmission  after  the  F-centers  are  formed  and 
relaxed,  it  is  found  that  there  is  no  change  in  the  two  photon  absorption 
coefficient.  In  order  to  determine  any  possible  absorption  other  than  two- 
photon  absorption  within  the  laser  pulse,  a  better  detection  system  to  reduce 
the  scattering  of  data  points  is  needed. 

The  polarization  dependence  in  KC1  has  been  measured.  A  change 
of  50%  is  observed  by  changing  the  electric  field  from  [00 1]  to  [Oil] 
direction. 

Since  the  power  used  is  much  higher  than  the  critical  power  for  self- 
focusing  [15],  it  is  important  to  notice  any  change  in  beam  size  by  self- 
focusing.  The  beam  sizes  before  and  after  the  sample  are  compared,  and 
the  measurements  have  been  done  with  different  sample  thickness,  also  by 
observing  the  stable  color  track  formed  in  U-center  doped  KC1  [16]'  under 
a  microscope.  No  evidence  of  self -focusing  has  been  observed  so  far. 

It  is  possible  that  n^  is  changed  by  the  two-photon  absorbing  levels,  and 
self-focusing  tends  to  concentrate  the  power  while  two-photon  absorption 
tends  to  remove  the  power  at  the  peak. 

An  electronic  system  to  monitor  the  energy  of  laser  pulse  has  been 
constructed,  but  noise  from  the  spark  gap  in  pockel  cell  single  pulse  selector 
presents  a  problem. 

Numerical  data  for  the  two  photon  absorption  coefficients  in  several 
alkali  halides  are  listed  in  Table  1.  These  should  be  compared  with  earlier 
data  listed  in  Table  2  but  obtained  by  less  accurate  techniques  available  a 
decade  ago.  These  were  also  quoted  in  the  Eighth  Technical  Report,  Dec¬ 
ember  1976,  by  M.  Sparks  and  C.  J.  Duthler  of  Xonics,  Inc.  It  is  seen  that 
the  intrinsic  two  photon  absorption  coefficients  appear  to  be  about  one  order 
of  magnitude  smaller. 


TABLE 
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It  is  planned  to  extend  the  measurements  of  the  two  photon 
absorption  cross  sections  to  other  materials  used  in  UV  windows  and 
coatings.  Furthermore,  data  at  the  third  harmonic  of  Nd-Yag  at  353  nm 
will  be  taken. 
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C.  ABSOLUTE  DETERMINATION  OF  THE  TWO-PHOTON 

ABSORPTION  COEFFICIENT  RELATIVE  TO  THE  INVERSE 
RAMAN  CROSS  SECTION 

With  the  development  of  high-power  lasers,  it  has  become  possible 
to  investigate  many  nonlinear  (NL)  optical  processes.  In  particular,  the 
two-photon  absorption  (TPA)  process  by  which  two  photons  are  simul¬ 
taneously  absorbed  has  been  explored  extensively.  The  selection  rules 
for  the  TPA  transitions  are  in  general  different  from  those  of  the  one- 
photon  transitions,  and  therefore,  the  measurement  of  the  TPA  properties 
of  a  material  system  is  complementary  to  the  linear  measurements.  A 
lot  of  information  on  the  symmetry  of  the  material  eigenstates,  including 

states  which  are  unseen  in  linear  spectroscopy,  may  be  obtained  from  the 

* 

relative  values  of  the  different  components  of  the  TPA  tensor.  Inform¬ 
ation  about  the  magnitude  of  the  transition  matrix  elements,  the  density 
of  states  and  the  energy  eigenvalues  may  be  obtained  from  absolute 

measurements  of  the  TPA  coefficient  and  especially  from  its  frequency 
2 

dependence.  Unfortunately,  these  latter  measurements  are  known  to  be 

difficult  because  of  the  strong  dependence  of  the  TPA  process  on  the 

parameters  of  the  lasers  used.  This  fact  is  reflected  by  the  large 

2 

scattering  of  the  reported  results  in  the  literature. 

New  techniques  for  absolute  measurement  of  the  TPA  coefficient 

3-9 

have  been  recently  reported.  In  the  first  kind  of  technique,  the  value 

obtained  for  the  TPA  coefficient  is  based  on  a  good  knowledge  of  the  laser's 
3  4 

parameters.  *  These  measurements  can  be  performed  only  with  those 
particular  lasers  which  may  be  quite  well  characterized,  and  therefore 
the  spectral  region  of  the  measurement  is  very  limited.  The  second  kind 
of  technique,  which  is  independent  of  the  laser  parameters,  is  based  on 
absolute  calibration  of  the  TPA  coefficient  by  other  known  nonlinear  (NL) 
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cross  sections.  This  type  of  measurement  may  be  carried  out  at  different 
laser  frequencies.  This  was  done  in  Ref.  5  using  the  second  order  NL 
susceptibility  of  quartz  as  the  NL  reference.  A  disadvantage  of  this 
specific  calibration  is  the  assumption  that  the  quantum  efficiency  for  the 
induced  fluorescence  by  the  absorption  of  one  photon  at  2 u  is  identical 
with  the  efficiency  associated  with  TPA  of  two  photons  at  u.  In  Refs.  6-9 
the  calibration  of  the  TPA  cross-section  is  done  in  relation  to  well  known 
Raman  scattering  cross  sections.  In  contrast  with  all  the  other  classical 
TPA  techniques ,  this  measurement,  which  utilizes  a  NL  parametric  mixing 
technique,  is  not  associated  with  real  material  transitions  and  therefore 

7 

stepwise  multiphoton  transitions  do  not  interfere  in  the  measurement 
and  the  value  for  the  TPA  coefficient  may  be  unambiguously  obtained. 

Another  simple  normalization  technique  has  been  proposed  recently*** 
which  is  very  useful  for  obtaining  relative  values  of  TPA  coefficients  of 
different  materials  independent  of  the  properties  of  the  laser  used.  **  This 
method  may  give  absolute  TPA  coefficients  by  using  reference  materials 
whose  TPA  coefficients  have  been  measured  by,  for  example,  the  method 
of  Ref.  7. 

The  purpose  of  this  paper  is  to  report  on  a  new  method  for  the  direct 

calibration  of  the  TPA  coefficient  relative  to  well  known  inverse  Raman 
8  12  13 

cross  sections.  ’  ’  By  this  method  we  calibrated  the  TPA  of  CdS, 

GaP,  SrTiOj  and  TiC^  at  31660  cm”*  relative  to  the  Raman  cross 
section  of  the  strong  2855  cm"*  CH  mode  of  liquid  cyclohexane.  Moreover, 
the  dependence  of  the  measurement  on  the  laser  parameters  could  be 
eliminated  in  a  two-channel  normalization  experiment  in  ways  similar  to  those 
used  in  Ref.  10,  and  thus  very  accurate  and  reliable  TPA  coefficients  were 
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obtained.  This  general  method  may  also  be  used  for  the  calibration  of 
absolute  Raman  cross  sections. 

Theory 

When  two  laser  beams,  a  strong  beam  at  u^  and  a  weaker  probe 

beam  at  Up  interact  in  a  NL  material  system,  the  strong  beam  induces 

a  change  in  the  complex  refractive  index  at  the  probe  frequency.  The 

real  part  of  this  change  is  proportional  to  nv  ,  the  self-focussing 

coefficient,  ^  while  the  imaginary  part  of  this  change  is  responsible 

14 

for  NL  loss  or  gain  at  Up.  The  net  gain  or  loss  at  Up  is  proportional 
to  the  imaginary  part  of  the  NL  third  order  susceptibility  x^^("wp»up> 
u l»  -  UjJ.  Neglecting  higher  order  processes,  the  probe  intensity  may  be 

4 

described  by: 

It  *  -  Vp  -  p'l'p  <u> 

where 

P  =  32  IT2  up/(c2npn;L)  Im  x(3Vup,  wp,uL, -(JL)  ,  (lb) 

1.  and  n.  {i  =  P,  L)  are  the  laser  intensities  and  the  refractive  indices, 
at  Up  and  u^,  respectively,  Up  is  the  linear  absorption  coefficient  at 
Up,  and  c  is  the  speed  of  light. 

It  is  clear  from  Eqs.  (1)  that  the  NL  process  is  enhanced  whenever 

X^  reaches  a  resonance.  ^  We  are  dealing  here  with  systems  which  are 

transparent  at  both  u^  and  Up,  and  therefore  only  the  resonances  of 
(3) 

X  which  are  associated  with  two-photon  transitions  (i.  e. ,  TPA  or  Raman) 
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of  the  material  system  will  be  taken  into  account.  For  a  system  occupying 

the  ground  state,  loss  at  Up  (or  a  positive  contribution  to  Imx^) 

occurs  in  the  TPA  case,  i.  e.,  when  up  +  u^  is  close  to  two-photon 

transitions  of  the  material.  Raman  transitions,  however,  may  contribute 

to  both  gain  and  loss  at  up.  When  u^  -  Up  is  close  to  a  Raman 

(3) 

transition,  and  >■  u p,  the  Raman  contribution  to  Imy  is  negative 

and  is  associated  with  gain.  This  gain  is  responsible  for  the  strong 

1 7 

amplification  in  the  stimulated  Raman  effect.  When  >  u^,  the 

(3) 

sign  of  the  Raman  contribution  to  Im*  is  positive  and  the  Raman 

transition  is  associated  with  absorption  at  Up.  The  latter  process,  which 

13  8  16 

is  usually  called  the  inverse  Raman  effect,  ’  ’  is  used  in  this  work 
for  the  absolute  calibration  of  the  TPA  coefficients.  The  cross  section  of 
the  inverse  Raman  effect  is  a  measurable  quantity  since  it  is  directly 
proportional  to  the  spontaneous  Raman  cross  section. 

We  will  now  explain  the  details  of  the  calibration  procedure.  We 
first  solve  Eq.  (1)  for  a  slab  of  NL  material  of  thickness  d  and  get  for 
the  transmitted  energy: 

Ip(x,  y,  z  =  d,  t)  =  (l-Rp)2Ip(x,  y,  z=  0,  t-np  d/c)X 

exp  -  {  ap  d  +  P^fljjx,  y,  z  =  0,  t-nL  d/c))  (2) 

Peff  =  P/qL  ’ 

where  p  is  the  TPA  coefficient  or  the  inverse  Raman  absorption 
coefficient  and  a.,  R^  and  n.  (i=P,  L)  are  the  one-photon  absorption 
coefficients,  the  reflectance  and  the  refractive  indices,  respectively.  In 
the  solution  (2)  we  assume  that  only  linear  attenuation,  and  no  back  reflection) 


affects  the  strong  beam  at  The  total  transmitted  energy  of  the 

probe  pulse,  E^,,  is  given  by  the  integration  of  Ip  over  the  pulse 
duration  and  the  pulse  cross  section.  Neglecting  the  dispersion  of  the 
refractive  index,  which  is  a  good  assumption  for  long  pulses  and  thin 
samples,  we  get  for  E^,: 

et  =  exp  { -apd}  (l-Rp)2jJdT  dx  dy  Ip(x,  y,  z  =  0,  t)  exp  -  {peffIL(x,  y,  z  =  0, 


This  expression  clearly  demonstrates  the  importance  of  the  spatial 
and  temporal  overlap  of  the  laser  pulses  in  the  measurement  of  p. 

Similarly  to  Ref.  10,  we  define  p(I^),  the  ratio  between  the  trans¬ 
mitted  energies  of  the  probe  beams  through  two  different  samples  A  and  B 
for  identical  laser  pulses: 


e*P(-a£  dA]  II I-d  exp  f  IT  }  dT  dx  dy 


exp{-apd  3/  /  exP  C-  P  eff  dT  dx  dy 


where  K  is  a  constant  independent  of  the  laser  intensity.  Note  that  p  is 
a  complicated  function  of  1^  and  the  beam  overlap.  If  the  strong  laser 
beam  is  attenuated  in  sample  A  by  a  filter  of  transmission  F,  the  factor 
1^  in  the  numerator  of  Eq.  (4)  should  be  replaced  by  F*  I^«  It  is  clear 
that  when  F  obeys  the  relation 

P  =  ?f>eAf£  •  <5> 


the  two  integrands  in  Eq.  (4)  are  identical  and  the  ratio  p  is  independent 
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of  1^  and  the  characteristics  of  the  beam  overlap.  Therefore,  since  F 

can  be  found  experimentally  by  checking  the  dependence  of  p  on  I^»  the 
B  A* 

ratio  Pejff/Peff  can  evaluated. 


Experiment  and  Results 

The  TPA  measurements  reported  here  were  calibrated  using  the 

2855  cm~*  mode  of  cyclohexane.  In  Fig.  1  we  illustrate  our  measurement 

of  the  inverse  Raman  absorption  cross  section  of  cyclohexane  in  the 

region  3000  >'  <jp  -  >2800  cm  \  This  measurement  was  done  by 

measuring  the  attenuation  of  the  variable  frequency  dye  laser  probe 
3  18  19 

beam  *  ’  7  for  a  normalized  intensity  of  the  strong  Ruby  laser  beam. 

With  the  narrow  linewidth  lasers  used,  the  measured  spectrum  of  Pg££ 

fits  quite  well  a  measurement  of  the  spontaneous  Raman  scattering  which 

is  also  presented  in  Fig.  1.  The  coincidence  between  the  two  spectra  is 

a  good  check  for  the  reliability  of  the  experimental  system. 

The  absolute  Raman  cross  section  of  the  2855  cm  *  mode  of  cyclo- 
.  o 

hexane  measured  in  the  blue  =  4880  A)  and  reported  by  Colles  and 

Griffiths^  is  used  for  calculating  x^(~Wp  >  wp»wL’  ~WL^  us*n§  t^le 
14  ^ 

relationship 


x(3,(- 


"WL) 


do 

d£2 


(6) 


(3) 

No  correction  for  the  small  dispersion  of  ^  was  introduced  in  the 
calculation  of  the  inverse  Raman  absorption  coefficient,  (3,  from 
The  cross  section  measurements  in  Ref.  20  are  based  on  the  Raman 
cross  section  of  the  992  cm  ^  mode  of  benzenereported  by  Skinner  and 
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Nilsen.  The  latter  cross  section  is  slightly  larger  than  the  one 

22 

reported  by  Kato  and  Takuma.  Therefore,  the  value  of  (3  for  cyclo¬ 
hexane  given  here  (p  =0.0019  cm/MW)  may  be  somewhat  high.  Con¬ 
sidering  also  experimental  reproducibility,  an  uncertainty  of  about 
±20%  for  the  cyclohexane  inverse  Raman  absorption  coefficient,  p,  seems 
to  be  reasonable. 

The  TPA  measurements  were  done  using  a  two-channel  normalization 

10  3  18  19 

technique  and  the  conventional  TPA  technique.  ’  ’  The  experimental 

setup  for  the  two-channel  normalization  measurements  is  shown  in  Fig.  2. 

A  multimode  Ruby  laser  is  used  as  a  source  for  the  strong  beam  at 

and  for  pumping  (after  doubling  the  frequency)  a  visible  dye  laser  (wp). 

An  intracavity  telescope  in  the  dye  laser  is  used  for  narrowing  the  output 
P  23 

below  0.  5  A.  The  beams  of  the  Ruby  and  the  dye  laser  are  combined 

by  a  dichroic-mirror  (DM).  A  glan  prism  (GS)  splits  the  colinear  beams 

into  the  cyclohexane  cell  (SA )  and  the  examined  sample  (S^).  After 

filtration  by  the  cutoff  filters  (F^),  the  two  probe  beams  are  detected  by 

two  HP  5082-4220  diodes.  A  plate  of  ground  glass  (G)  is  inserted  in 

front  of  each  diode  in  order  to  minimize  the  effect  of  the  nonuniformity  of 

the  laser  beam  and  the  diode's  active  area.  The  signals  are  amplified, 

9 

digitized  and  punched  on  a  tape  which  is  later  analyzed  by  a  computer. 

The  relative  intensity  of  the  Ruby  laser  beam  in  the  samples  (i.  e.  ,  the 
value  of  F)  may  be  varied  by  rotating  the  double  Fresnel  prism  (DFP)  used 
as  a  half-wave  plate.  For  a  given  value  of  F  the  ratio  p  (IjJ  is  measured 
for  different  Ruby  intensities  which  are  obtained  by  inserting  filters  (VA) 
in  the  Ruby  beam.  An  iris  in  front  of  the  splitting  prism  GS  is  a  critical 
component  in  the  TPA  measurement.  By  proper  setting  of  this  iris  the 
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overlap  integral  [Eq.  (3)]  is  maximized. 

Fig.  3a  shows  our  results  for  p  for  CdS/cyclohexane  and  for 

CdS/TiC^*  It  is  seen  from  the  figure  that  the  appropriate  values  of  F 

for  which  p  is  independent  of  1^  are  2.  2  <  F  <  2.  56  and  Fa  1.4 

for  CdS/C^H^  and  CdS/TiO^,  respectively.  Using  Eq.  (5)  and  the 

linear  data  of  the  samples  given  in  Table  I,  we  obtain  the  values 

^CdS  =  i  •  006  cm/MW  and  P^dS^TiO  2.  4.  In  a  similar  way 

2 

we  normalized  (3  of  SrTiO^  vs.  cyclohexane,  PgrTiO  =  .003  ±.  001  cm/MW. 

Using  the  conventional  TPA  technique  we  measured  the  TPA  coeffi¬ 
cients  of  CdS,  TiO^,  and  GaP  relative  to  the  inverse  Raman  cross  section 
of  cyclohexane.  The  TPA  coefficients  were  obtained  by  comparing  the 
transmission  of  the  samples  at  Up  for  different  values  of  1^.  In  order 
to  minimize  the  drift  in  the  experimental  conditions,  the  samples  were 
mounted  on  a  special  holder  which  allowed  us  to  switch  samples  rapidly. 

The  results  of  this  measurement  are  shown  in  Fig.  3b  and  are  tabulated 
in  the  fifth  column  of  Table  I.  For  comparison  the  results  of  the  two- 
channel  normalization  measurements  are  listed  in  the  fourth  column  of 
this  table.  The  error  bars  in  Fig.  3  represent  the  standard  deviation  of 
p  and  Pe££*  The  uncertainty  in  the  data  points  for  cyclohexane  and  TiOz 
shown  in  Fig.  3b  are  relatively  small  because  of  the  large  number  of  laser 
shots  (~50)  averaged  at  each  of  these  data  points. 

The  value  P^^  =  .071  ±.015  cm/MW  and  the  ratio  PcdS^TiO  = 

2 

2.  6  ±  .  6  obtained  by  the  conventional  TPA  technique  fits  very  well  the 
results  of  the  two-channel  normalization  measurement.  The  uncertainty 
in  the  former  measurement  is,  of  course,  larger  than  that  of  the  two- 
channel  experiment. 


( Arb.  Units ) 


FIGURE 


a)  p,  the  ratio  between  the  transmitted  energies  in  the  samples 
CdS/cyclohexane  and  CdS/TiO^  vs.  Ruby  laser  intensity.  Each 
value  of  F  corresponds  to  a  setting  angle  a  of  the  DFP. 

b)  vs.  Ruby  intensity,  as  we  measured  by  the  conventional 
TPA  technique.  Each  data  point  represents  the  average  of 

^  10  laser  shots  except  the  TiO-  and  cyclohexane  with  ~50 
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Our  measured  value  for  0  of  TiO->  at  u>.  +  oj—  =  31660  cm 

ZZ  C  la  Xr 

is  about  7  times  larger  than  the  value  (3  =  .  0034  cm/MW  obtained  by 

z  z 

24 

Waff  and  Park  (who  used  a  Nd-glass  laser  and  flash  lamp)  at  the  same 

two-photon  energy,  and  about  3  times  larger  than  {3  =  .  0065  cm/MW, 
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obtained  by  Penzkofer  and  Falkenstein  at  2u  =  31950  cm  using  a  pico¬ 
second  laser  source.  (6  and  8  represent  the  TPA  coefficient  for 
linearly  polarized  photons,  polarized  along  the  x  axis  and  the  z  axis  of 
the  crystal,  respectively).  No  TPA  data  at  +  Up  ~  32000  cm  *  for 
the  other  crystals  investigated  in  this  work  has  been  published  in  the 
literature. 

A  comparison  of  our  TiO^,  results  with  those  of  Refs.  24  and  25 

should  be  carefully  made.  We  note  that  the  most  fundamental  physical 

(3) 

quantity  to  be  compared  is  x  and  not  P>  which  includes  a  trivial 
frequency  factor  [see  Eqs.  ( 1 ) J.  The  large  deviation  of  the  result  of 
Ref.  24  probably  originates  from  the  severe  complications  which  arose 
in  their  experiment  in  the  determination  of  the  overlap  integral  [Eq.  (3)]. 
The  deviation  of  the  results  of  Ref.  25,  however,  cannot  be  explained  by 
the  different  laser  frequencies  used,  by  the  small  difference  between  the 
two-photon  energies  or  by  the  different  electric  field  polarizations  used, 
which  according  to  Ref.  24  result  in  the  ratio  (3  /|3  =  1.25.  Dispersion 

ZZ  XX 

in  the  TPA  coefficient  which  originates  from  the  one -photon  denominators 
(3) 

in  x  is  also  unlikely  as  an  explanation  for  this  deviation  since  the 

measurements  were  done  with  laser  frequencies  far  enough  from  the  band 

gap  (3.  75  eV)^  to  insure  no  resonant  behavior.  ^  Also,  no  [x^]^ 

(3)  12  9 

contribution  to  x  >  ’  a  contribution  which  depends  on  the  beam 

polarization  and  propagation  vector,  is  expected  in  TiO^  which  exhibits 
inversion  symmetry.  Two  steps  sequential  one-photon  absorption 


-30- 


contributions  to  p  and  absorption  caused  by  two-photon  generated  holes 
may,  however,  influence  differently  the  measurements  in  the  nanosecond 
and  picosecond  regimes.  ’  *  We  conclude,  therefore,  that  the  latter 
effects  and  the  possibility  that  the  laser  pulses  in  Ref.  25  were  not  well 
characterized  may  explain  the  deviation  of  their  TiO^  TPA  results  from 
ours. 


Discussion  and  Conclusions 

The  TPA  coefficient  and  the  Raman  cross  section  are  important 

material  parameters.  Accurate  values  of  TPA  coefficients  and  especially 

their  spectral  dependence  will  be  useful  for  stimulating  further  theoretical 

calculations  of  two-photon  phenomena.  In  particular,  existing  theoretical 

27 

models  of  matter  such  as  the  "empirical  pseudopotential  method"  may 
be  checked  by  comparing  the  theoretical  p  based  on  the  specific  model 
with  the  experimental  one.  We  note  that  one-photon  spectroscopy  above 
the  material  band  gap  show's  surface  effects  because  of  the  small  pene¬ 
tration  depth.  The  TPA  measurement,  however,  reflects  in  the  same 
spectral  region  only  the  bulk  properties  of  the  material.  Therefore,  a 
comparison  of  the  TPA  experimental  results  with  theoretical  models  is 

more  direct  than  in  the  one -photon  case.  Experimentally,  accurate  TPA 

25 

coefficients  may  be  used  for  the  characterization  of  pulsed  laser  beams. 
This  may  be  done  by  checking  the  interaction  of  the  strong  beam  with  a 
probe  beam  whose  parameters  are  well  known  (cw  dye  laser  beam,  for 
example)  in  a  NL  absorber.  The  measured  value  of  E^,  [see  Eq.  (3)] 
may  give  important  information  on  the  spatial  and  temporal  properties  of 
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the  strong  beam.  In  particular,  absolute  peak  power  values,  which  are 

very  important  in  multiphoton  processes,  may  be  obtained  by  this  method. 

Accurate  Raman  cross  sections  may  be  used  for  measuring  average 

concentrations  of  molecular  species  over  long  atmospheric  paths  using 

19 

stimulated  Raman  gain  or  inverse  Raman  loss  experiments. 

The  calibration  method  described  here  is  very  useful  for  reliable 
measurements  of  TPA  coefficients.  Raman  cross  sections  may  also  be 
calibrated  by  the  same  method.  In  this  case,  inverse  Raman  absorption 
as  well  as  stimulated  Raman  gain  may  be  utilized  in  the  measurement. 

The  use  of  a  probe  beam  in  the  absorption  measurement  has  an  advantage 
over  the  standard  one-beam  TPA  methods.  In  the  one-beam  case,  any 
small  loss  in  the  laser  beam  caused  by  other  NIL  processes  may  strongly 
influence  the  measured  cross  section.  Such  perturbations  which  influence 
the  strong  laser  beam  have  only  a  minor  effect  on  the  probe  beam. 

In  conclusion,  we  have  described  a  technique  which  allows  simple 
measurement  of  TPA  coefficient  and  Raman  cross  section,  relative  to 
known  inverse  Raman  cross  sections.  The  use  of  the  two-channel 
normalization  experiment  enables  one  to  eliminate  the  dependence  of  the 
NL  absorption  on  the  laser  intensity,  thereby  removing  the  usual  problems 
associated  with  absolute  intensity  measurements.  Additionally,  this 
normalization  alleviates  sources  of  errors  due  to  irreproducible  laser 
pulse  spatial  and  temporal  profile,  allowing  the  measurement  to  be  made 
with  simple  uncontrolled  lasers. 
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D.  NONLINEAR  ABSORPTION  AND  REFLECTION  IN  InSb 

DURING  SHORT  COz -LASER  PULSES 

Nonlinear  interactions  of  light  in  semiconductors  has  been  studied 
intensively  over  the  last  fifteen  years  by  various  authors  [1-3].  In  the 
high  intensity  limit,  the  material  of  course  breaks  down  and  suffers 
irreversible  damage  [4],  In  the  low  intensity  limit,  a  perturbative  approach 
is  usually  taken  whereby  the  interaction  is  described  as  a  multiphoton 
process  [5],  It  is  in  the  region  just  below  the  material  breakdown  threshold 
that  is  of  fundamental  physical  interest  and  is  relatively  unexplored  experi¬ 
mentally.  In  this  regime,  electron-hole  generation  is  unavoidable  and  the 
interaction  of  light  with  the  material  is  a  mixture  of  intrinsic  material 
interaction  and  plasma  effects. 

In  this  section  a  study  of  the  behavior  of  InSb  under  interaction  with 
high  intensity  10.  6  radiation  is  described.  In  the  past,  experimental 
studies  of  plasma  behavior  were  done  with  picosecond  1.06  pulses  [6-8], 
This  is  because  of  the  ability  to  go  to  high  peak  power  without  damaging 
the  crystal  since  damage  thresholds  are  mostly  energy  thresholds.  We 
are  able  to  generate  50  psec.  C02  pulses  with  ~  0.  5  Mw  power  and 
therefore  can  extend  such  studies  to  small  band  gap  materials  [9].  In 
addition  to  being  a  pulse  shaper  and  limiter  for  high  power  lasers  [10], 

InSb  is  interesting  because  of  its  highly  nonparabolic  energy  bands  and 
associated  nonlinearities.  Its  behavior  at  high  intensities  has  been  studied 
theoretically  by  Tzoar  and  Gersten  [  1 1  ] .- 

Figure  1  shows  the  transmitted  energy  of  a  sample  of  InSb  doped 
16  -  3 

to  n  =  6.  6  X  10  cm  versus  incident  laser  intensity.  It  is  seen  that 
the  relative  transmittance  decreases  by  a  factor  of  four  at  high  intensities. 
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This  is  probably  due  to  the  generation  of  a  critically  dense  electron-hole 
plasma  rather  than  by  two-photon  absorption  alone.  To  confirm  this,  we 
made  a  calculation  similar  to  Bechtel  [12]  to  determine  the  two-photon 
absorption  coefficient.  From  our  experimental  data  for  the  transmittance 
a  value  of  9.  1  cm/Megawatt  was  derived,  which  is  a  factor  of  35  larger 
than  the  established  value  of  0.  23  cm/Mw  [10,  13].  The  reason  is  that 
carrier  effects  cannot  be  ignored  in  the  calculation.  In  view  of  this,  it  is 
more  advisable  to  look  at  the  reflectance. 

Figure  2  shows  the  reflectance  of  a  highly  doped  sample  (  n  =  1.  5 
18  -  3 

X  10  cm  )  versus  intensity.  A  brewster  angle  of  incidence  was  used. 
The  reflectivity  changes  by  more  than  a  factor  of  six.  What  we  actually 
measured  was  the  integrated  transient  reflectivity  R(t)  [14].  The 
increase  in  reflectivity  is  a  manifestation  of  the  changing  index  of  refrac¬ 
tion  due  to  formation  of  a  plasma.  The  electron  generation  mechanism 
is  by  multiphoton  injection  [13,  15], 


3n  _  o_I  +  3l^  n 
3t  tiu  2  ftu  ~  t 


yn(n+nQ) 


where  y  is  the  radiative  recombination  rate  and  t  is  the  nonradiative 
recombination  time  which  is  mostly  due  to  Auger  recombination,  a  is 
single  photon  excitation  from  intermediate  levels  and  8  is  the  TPA 
coefficient.  The  single  photon  term  has  been  measured  by  Fossum  [13] 
at  77°K.  With  their  number  for  a  and  at  our  intensity,  the  first  term 
is  four  orders  of  magnitude  smaller  than  the  second  term.  It  is  expected 
to  be  even  smaller  at  high  temperature  since  the  intermediate  states  are 
all  thermally  ionized  to  begin  with.  The  estimated  recombination  lifetime 
is  500  nsec,  and  the  nonradiative  recombination  time  is  2  3  nsec.  [10»  13]. 


REFLECTED  SIGNAL  (v) 


INCIDENT  SIGNAL  (xl.2  xl09W/cmz) 


FIGURE  2. 


Reflected  energy  versus  incident  intensity  for 
1]  =  1.  5X  10 cm"^  sample.  Black  dots:  data 
points  for  reflection  from  a  polished  aluminum 
mirror.  Open  dots:  data  points  for  sample. 
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Therefore,  recombination  is  largely  nonradiative.  Since  both  times  are 
much  larger  than  the  pulsewidth,  they  can  be  ignored,  and  the  equation 
describing  plasma  generation  is  simply 

3n  _  o  I2 
St  P  2  hu 

For  small  reflectivities  in  the  region  of  the  brewster  minimum,  it  can 
be  shown  that 


However,  the  experimental  curve  in  Figure  2  shows  that  R  is  approximately 

proportional  to  the  intensity,  implying  that  n  grows  only  as  the  square 

root  of  the  intensity.  This  is  ah  indication  that  our  simple  model  of  two 

photon  ionization  is  not  valid  at  such  high  intensities.  A  straightforward 

Keldysh  type  calculation  [l,  16]  shows  that  the  ionization  process  changes 

from  a  multiphoton  character  to  a  tunneling  character  at  around  1 ^ 

9  2 

=  5.  5  X  10  W/cm  for  our  sample.  In  the  tunneling  regime,  the  rate 
of  carrier  generation  follows  a  Keldysh  type  of  equation 

«  I2(Icr/I)3/4exp{-3.9(Icr/I)1/2(l-Icr/8I)} 

This  expression  is  in  rough  agreement  with  the  experimental  trend.  A 
quantitative  calculation  should  take  into  account  the  temporal  variation 
of  I  during  the  pulse.  The  accuracy  of  the  measurement  must  also  be 
further  improved. 


-40- 


References 


1.  L.  V.  Keldysh,  Zh.  Eksp.  Fiz.  47,  1945  (1964)  [Sov.  Phy.  Jet-P 
20,  1307  (1964)]. 

2.  Man  Mohan  and  S.  N.  Hague,  Phys.  Rev.  A9  ,  631  (1974); 

G. M.  Aratyunyan  and  S. M.  Shahinyan,  Phys.  Stat.  Sol.  (b)  77, 

K171  (1976). 

3.  Howard  R.  Reiss,  Phys.  Rev.  A 1,  803  (1970). 

4.  W.  Lee  Smith,  J.H.  Bechtel  and  N.  Bloembergen,  TR  No.  665, 
Division  of  Engineering  and  Applied  Physics,  Harvard  University. 

5.  R.  Brauns tein  and  N.  Ockman,  Phys.  Rev.  134,  A499  (1964). 

6.  C.  V.  Shank  and  D.  H.  Auston,  Phys.  Rev.  Lett.  34,  479  (1975). 

7.  John  F.  Reintjes  and  James  C.  McGroddy,  Phys.  Rev.  Lett.  30, 

901  (1973). 

8.  D.  J.  Kennedy  et  al.  ,  Phys.  Rev.  Lett.  32,  419  (1974). 

9.  H.S.  Kwok  and  E.  Yablonovitch,  Appl.  Phys.  Lett.  30,  158  (1977). 

10.  A.F.  Gibson  et  al.  ,  J.  Phys.  CC:  Solid  State  Phys.  9.,  3,  259  (1976). 

11.  N.  Tzoar  and  J.I.  Gersten,  Phys.  Rev.  B12,  1132  (1975). 

12.  J.H.  Bechtel  and  W.  L.  Smith,  Phys.  Rev.  13,  3515  (1976). 

13.  H.  J.  Fossum  and  D.  B.  Chang,  Phys.  Rev.  8,  2842  (1973); 

H.  J.  Fossum,  W.S.  Chen,  and  B.  Ancker-Johnson,  Phys.  Rev. 

138,  2857  (1973). 

14.  S.A.  Jamison,  A.V.  Nurmikko,  and  H.  J.  Gerritsen,  Appl.  Phys. 
Lett.  29,  64  (1976). 

15.  R.  E.  Slusher,  W.  Giriat  and  S.  R.  J.  Brueck,  Phys.  Rev.  183, 

758  (1969). 

16.  S.S.  Mitra,  L.M.  Narducci,  R.A.  Shatas,  Y.  F.  Tsar  ind  A. 
Vaidyanathan,  Appl.  Opt.  _f4,  3038  (1975). 


